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Abstract 

Effect of ionic solute on a near-critical binary aqueous mixture confined between charged walls 
with different adsorption preferences is considered within a simple density functional theory. For 
the near-critical system containing small amount of ions a Landau-type functional is derived based 
on the assumption that the correlation, ^, and the Debye screening length, k~^, are both much 
larger than the molecular size. The corresponding approximate Euler-Lagrange equations are 
solved analytically for ions insoluble in the organic solvent. Nontrivial concentration profile of the 
solvent is found near the charged hydrophobic wall as a result of the competition between the 
short-range attraction of the organic solvent and the electrostatic attraction of the hydrated ions. 
Excess of water may be present near the hydrophobic surface for some range of the surface charge 
and i^K. As a result, the effective potential between the hydrophilic and the hydrophobic surface 
can be repulsive far from the critical point, then attractive and again repulsive when the critical 
temperature is approached, in agreement with the recent experiment [Nellen at.al.. Soft Matter 7, 
5360 (2011)]. 
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Near-critical binary mixture confined in a slit induces effective attraction or repulsion 
between the confining walls if adsorption preferences of the two walls are the same or oppo- 



site respectively [lll3|. The range of this so called thermodynamic Casimir force is of order 
of the bulk correlation length ^. Parallel walls covered by electric charges of the same sign 
repel each other. One could thus expect stronger repulsion between likely charged walls with 
opposite adsorption preferences confining near-critical binary mixture. In striking contrast 



to the above expectation, in the recent experiments j4| strong attraction was observed be- 
tween a charged hydrophobic colloid particle and a charged hydrophilic substrate for some 
range of temperatures and concentrations of a hydrophilic salt added to the solution. Ef- 
fective interactions between the colloid particles separated by distances much smaller than 
their radii are similar to the interactions between planar surfaces. Possibility of changing 
these interactions from attraction to repulsion by minute changes of temperature or salinity 
opens possibilities for designing and controlling reversible structural changes, in particular 
aggregation or adsorption. It is thus important to understand the mutual influence of the 
critical adsorption and the distribution of ions that leads to the attraction between the walls 
instead of the expected repulsion. We address this issue in this communication. 

We consider a water - organic liquid mixture containing hydrophilic ions in a slit with 
selective, charged walls of the area A — )• oo, separated by the distance L ^ 1 (Fig.l). We 
choose the average diameter of the molecules, a = 1, as the length unit and all the corre- 
sponding functions are dimensionless. The grand thermodynamic potential of the system 
can be written in terms of the local dimensionless densities Pi{r), where i = 1,2,3,4 for 
water, oil, + and — ion respectively, in the form 

n = -pAL + AUex (1) 
= iJ^d^ X - r')g^j{r - r')p,(r') + drp,{r){Vf{r) - fi,) + U^i - TS 

where p is the bulk pressure, the integration is over the system volume V = AL, periodic 
boundary conditions are imposed in the directions parallel to the walls, and S, Uei, T and 
/ij are entropy, electrostatic energy, temperature and chemical potential of the z-th specie 
respectively. Vij and gij are the van der Waals (vdW) interactions and the pair correla- 
tions between the corresponding components respectively, and the summation convention 
for repeated indices is assumed in the whole communication. V^'*(r) is the sum of the direct 
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wall- fluid potentials acting on the component i. Finally, 



= (7o + 7l) + ^(^) (2) 

is the excess grand potential per surface area, 7„ is the surface tension at the n-th wall 
{n = 0,L), and the effective potential "^{L) is the subject of our study. Because of the 
translational symmetry in the parallel directions, the densities depend only on the distance 
from the left wall, z. We make the standard approximation 

S/A = -kby2 dzpi{z)\npi{z), (3) 



and the standard assumption 







U4(j)]/A = r dz[-^(vV')' + e(j)ij] + ea(0)^(0) + ea(L)V'(L), (4) 
Jo on 

where the electrostatic potential ip satisfies the Poisson equation, 

e d'^ip{z 



+ e^z) = 0, (5) 



e is the elementary charge, e is the dielectric constant of the solvent, cr{n) is the dimensionless 
surface charge density at the n-th wall, and 

(j){z) = p3{z) - p^iz) (6) 

is the dimensionless charge density. Compressibility of the liquid can be neglected, and we 
assume J2i=i Pi = 1- We choose 0, the solvent concentration s = pi — p2, and the density of 
ions Pc = P3 + P4: as the three independent variables. Bulk equilibrium densities for given T 
and Pi correspond to the minimum of —pAL, and are denoted by s and pc- In equilibrium 
(f){z) and the deviations from the bulk values, 

'dliz) = pi{z) - p2iz) - S, ^2{z) = p3{z) + piiz) - pc (7) 

correspond to the minimum of C0ex['^i,i^2, <P] = {^['^i + s,'d2 + Pc,(p] — ^[s, Pc,0])/A with s 
and Pc fixed. We choose for s and pc the values corresponding to the critical point. 

Common salts are soluble in water and insoluble in organic solvents. We thus assume 
that the difference in the chemical nature of the anion and the cation is negligible, and 
postulate the same vdW interactions, V^^s = Vi^^. From Eq.([T]) it easily follows that the vdW 
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FIG. 1: Model system consisting of water, organic liquid (for example lutidine) and ions between 
negatively charged hydrophilic (dark, blue) and hydrophobic (light, red) walls. 

contribution to the internal energy expressed in terms of the new variables is independent of 
when Vi^3 = V^^4 [6|. Because Uei is independent of -di (see dl])), in this approximation the 
vdW and the electrostatic contributions to the internal energy are decoupled. The coupling 
is present in the entropic part. The excess entropy, Sex['&i,^2, <P] = {Sldi + s,^2 + Pc,<P] — 
S[s, pc, 0])/A, can be Taylor expanded in terms of 'i9j(;z) and (piz). For a near-critical mixture 
with small amount of ions the expansion can be truncated, because i^ilz) and (f){z) are small 
(except from microscopic distances from the surfaces). Using Eq.(l3]) one can verify that the 
excess entropy contains no terms proportional to and the lowest-order mixed term 

is (j)'^'d2- Thus, the excess grand potential can be split in two leading-order terms and the 
correction AC 

From the minimum condition for Uex it follows that the linear terms vanish, and the dominant 
terms in Eq.([H]) are quadratic in the fields '(9j(-2) and (f){z). The second term on the RHS of 
Eq.dH]) has the form 

Cdh[<P] = UM/A+^ f'dz{<p\z) + 0{<p')). (9) 

2pc Jo 

Eqs. dH]), dl]) and ([5]) agree with the Debye-Huckel (DH) theory for the excess grand potential 
of ions in a homogeneous solvent confined in a slit with charged walls. The first term in 
Eq.dS]) is equal to the excess grand potential per unit area for one kind of neutral solute in a 
two-component solvent, where the excess concentration of the solvent and the excess solute 
density are denoted by "i^i and -(92 respectively and the total density is fixed. This is because 
we assumed no difference between the vdW interactions of the anion and the cation - when 
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uncharged, they represent the same species in this theory. Close to the critical temperature 
Tc the fields ^i{z) vary on the length scale ,^ oc| (T — Tc)/Tc l^'^^ 1 with z/ ^ 0.63, and the 
standard coarse-graining procedures leading to the Landau functional can be applied |8|. 
Our coarse-graining of the first term in Eq.([T]) (expressed in terms of the new variables) is 
based on the Taylor expansion of "^i^z') about z' = z. The excess grand potential is expressed 
in terms of the fields and their derivatives, and in terms of the appropriate moments of 
the vdW interaction potentials, 

+mMM _ ,,(o),.(o) + mMM _ 



where 



= - 7° - ^ ^^ex [^1,^2,0] 



(11) 



J°. = / drJijir) and Jij = | / drJij{r)r'^. —Jij{r) represents the vdW interactions for -di and 
'dj, and can be obtained from the vdW contribution to Eq.([l]) with the densities expressed in 
terms of the new variables. We shall assume that the interaction ranges Qj defined by = 
6Jij/J^j are all Qj ~ 1, and characterize the system by three interaction parameters, Jij = 
Jij/6 (for the length unit a = 1). Finally, hi{n) is the surface field describing interactions 
with the n-th wall. The remaining surface terms result from the compensation for the 
interactions with the missing fluid neighbors at the wall. These interactions are included in 
the bulk term, but should be absent if the wall is present. In Ref. 6| the same functional 
was obtained from a lattice model for the four-component mixture. When the mixture 
phase separates, both the solvent concentration and the density of solute are different in the 
coexisting phases, because the solute is soluble only in water. Thus, Cc must depend on 
both, -^i and ^2- The critical order parameter is the eigenvector of C^j corresponding to the 
eigenvalue vanishing for T = Tc. 

In Eqs. lfTOj) and ([9]) the terms proportional to ksT represent the leading-order contribu- 
tions to the excess entropy per surface area. The next-to-leading order contribution has the 
form 



AC[d2,(P] = -ksT ["^dz 
Jo 



2p: 



o2 



(12) 
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and results from the fact that there are more ways of introducing a local difference in the 
concentrations of the anions and the cations, (plz), in the regions where there is more ions 
{'&2{z) > 0) than in the regions where there is less ions {'&2{z) < 0). 

The Euler-Lagrange (EL) equations for the functional - f[T^ . with the higher order 
terms in ilTUI) . and (IT^ neglected, take the forms 



dz^ 



M,,d,{z)+d4\z) (13) 



^ = "^^'^ + — d? — • 

The charge neutrality condition, dz(j){z) + ao + aL = 0, is imposed on 0, and the boundary 
conditions for i)i are 

^^U=o-MO) = Hm (15) 

In the above Mij = {J~^)ikC^j, where {J^^)ik is the (z, A;)-th element of the matrix inverse 
to the matrix Jij [6]. The remaining parameters are (rfi, (^2) = —^^(^{J~^)i,2, {J~^)2,2^, and 
Hi{n) = {J~^)ijhj{n). For a hydrophilic (hydrophobic) wall i^i < {Hi > 0). 

When A£ in Eq.([8]) is neglected, the Casimir and the electrostatic potentials are inde- 
pendent contributions to Wex, and the EL equations are linear and decoupled (the second 
terms on the RHS of Eqs. (fT3l) and (JT^ are absent). In a semi-infinite system the solutions 
of the linearized EL equation ( |T4l) and ( |T3l) are 



where = Q-Iq 



(f)^^\z) = —Kaexp{—Kz) (16) 

9|], and 

= Aexp{-z/0 + Qexp{-X2z), (17) 



where Ai and Cj depend linearly on Hi. The superscript (1) refers to the solutions of the 
linearized EL equations. In the critical region ^ ^ 00 and A2 ^ 1/^, therefore the second 
term on the RHS of Eq. fll7l) can be neglected. In the slit the equilibrium fields (p'^^^z) and 
'&'l^\z) contain also terms oc exp(— fi:(L — z)) and oc exp(— (L — z)/^) respectively (and the 
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amplitudes are modified). The effective potential "^{L) is obtained by subtracting the L- 
independent part from Uex calculated for the equilibrium profiles. Neglecting AC in ([8]) we 
obtain $(L) = dh{L) + ^c{L) with ^ dh{L) oc exp(-/tL) and ^c{L) oc exp(-L/^). 

The nonlinear terms in the EL equations f[T^ and f[T^ can be neglected when their 
magnitudes are much smaller than the magnitudes of the linear terms at the relevant length 
scales. The nonlinear contributions to Eqs. (fT3|) and can be estimated by examining 

z = K respectively. This is because the linear terms 
in Eqs. fll3l) and (IT^ decay on the length scales ^ and respectively. From Eqs. fll6l) and 
ffT7|) we obtain (p^'^^C) exp(— ^/t) and 'd^2\'^^^) ^ 6xp(— (^k)^^). Thus, the magnitudes of 
the correction terms depend crucially on the ratio between the correlation and the screening 
length, ^K. When — )■ oo, then 0(1) (0 and d^2\^^~^) = 0(1), therefore we may 
consider linearized f|T3|) . and treat f|T^ perturbatively. This case was considered in Ref. 6| 
for a semi-infinite system, and in Ref. (iQj for a slit. On the other hand, for — > we 
have 0(1) (0 = 0(1) and 

{k 1) —7- 0, therefore linearized Eq. ffMl) can be considered. As 
a consequence, 0^ in Eq. (IT^ can be approximated by 0(i)^. In this approximation Eq. (IT5]) 
takes the form of a linear inhomogeneous equation. We assume that this approximation is 
reasonable as long as .^/t < 1, and the magnitudes of cr^ and Ai are comparable. 

In the experiments showing unusual dependence of the effective potential '^{L) on T, 



and consequently on ^/t, the relevant lengths ratio was < 1 4", therefore in this work 
we assume = Eq. (fT3|) with approximated by 0*-i) can be easily solved analytically. 
The excess concentration of the solvent in the semi-infinite geometry takes the form 



d,{z) 



exp{-z/0 (18) 



where Bi = -^kBT, B* depends on the vdW interactions and 



h{K, y) = ^ f{y){2y - 1) - (1 + 2K)Ky 
y + k\ 



B* ^f{y) 



1+1 



(20) 



The excess solvent concentration at the distance z from the hydrophobic surface with weak 
and strong surface charge is shown in Fig. 2 for a few values of .^/t < 1 (note that in the 



figure captions the length unit a is re-introduced). In all the cases we observe excess of 
organic liquid close to the surface. In some cases, however, di{z) is non-monotonic and 
changes sign for zq ~ ^. Excess of water appears at the distances z > zq from the surface 
for all values of .^k < 1 in the case of strong surface charges. For weak surface charges 
excess of water appears only for < yo{<7); for > yo^cr) a monotonic decay of 
occurs, as in the uncharged system. Thus, the presence of the surface charge can change 
a (weakly) hydrophobic surface to an effectively hydrophilic one if we pay attention to the 
concentration of water at sufficiently large distances from the wall, 2 ~ ^. Chan ge o f the 
adsorption preferences by increased surface charge was observed experimentally 0, 11|. We 
emphasize that the change of the adsorption preference for small or moderate surface charges 
is present only sufficiently far from the critical point. 

The above properties can be understood by examining Eq. f|T8|) for the hydrophobic surface 
{Ai < 0). For simplicity we neglect ifii^^'^O — for ^ — )• 00 and k — )■ i.e. for 
T ^ Tc and pc — 0), and assume > 0.5. For > 0.5 the second term in Eq.l fTSi) decays 
faster, and at the length scale ^ the excess of water is found when the prefactor of the first 
term is positive. From Eqs.( |T8|) and ( |T9i) we can conclude that for T — )■ Tc (i.e. ^/t — )■ 00) 
the excess of water occurs (i.e. > 0) when Ai + ^-Bi > 0, which leads to the condition 

for the surface charge > 4|yli|/i?i. When cr^ < 4|Ai|/_Bi, excess of organic liquid occurs 
(i.e. < 0) for T — Tc. Since a'^Bif^^n) increases substantially when decreases (see 

Eq.(!T9l)). the prefactor of the first term in Eq. lfTSl) . Ai + a'^Bif{^K), changes sign for = yo- 
Thus, for 0"^ < 4|Ai|/i?i a crossover from the excess of the organic liquid for > (close 
to Tc) to the excess of water for < y^ (far from T^) occurs for sufficiently large distances 
from the hydrophobic surface, ~ ^. 

Physics behind such behavior is quite simple. The charged wall with no adsorption 
preference attracts ions. The ions insoluble in the organic liquid attract in turn water 
molecules to this wall. The excess number density of the hydrated ions (and thus the 
excess of water) appears in the layer of the thickness ~ {2k) and depends on the 

surface charge. The charge-neutral, hydrophobic surface attracts organic molecules. Excess 
of organic liquid is found in the layer of thickness ~ ^, and depends on the hydrophobicity of 
the surface. Competition between the excess of organic liquid and the excess of water near 
the surface which is both hydrophobic and charged depends on on the surface charge 
and on the hydrophobicity of the wall, and leads to the nontrivial concentration profiles. 
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FIG. 2: The excess solvent concentration •di{z) at the distance z from the charged hydrophobic 
surface for different values of (,k, with k = O.la^^, where a is the molecular diameter. The dimen- 
sionless surface charge density (Eq.Q) is (a) a = 0.048/a2 and (b) a = 0.092/0^, and in Eqs.([l8D 
-([20]) Ai = -0.22a-3, B* = OM{kBTca'^)-\ {J~^)i,2/B* = -12.52 and = 1.08 • lO^^a-s. For 
a = Inm (approximate size of the lutidine molecule) pc ~ 1.8 • lO^'^lmol/lit]. z and i^i are in a 
and units respectively and '&i{z) > for excess of water. 

The Casimir potential between the walls results from the change of the concentration 
near the first wall caused by the presence of the second wall. Let us consider vicinity of the 
hydrophilic wall when the weakly hydrophobic wall is present at the distance L ~ ^. The 
uncharged hydrophobic wall leads to depletion of water, but as discussed above and shown in 
Fig. 2, in the presence of the surface charge the hydrophilic ions can lead to the opposite effect. 
Thus, for the range of temperatures corresponding to the change of the adsorption preference 
of the weakly hydrophobic surface, the Casimir potential can be attractive. For not too 



9 




FIG. 3: The effective potential per unit surface area between the charged hydrophilic and hy- 
drophobic surfaces for the model system with J^^ = IksTca^ ^ — J22 ~ O-Olfc^Tca^, for 
K, = O.la^^ and different values of shown in the inset. The surface fields (Eq. psp l are 
Hi{{}) = -0.003a-3,i7i(L) = 0.002a-3, i72(0) = -H2{L) = -O.Sa-^, = 1.08 • IQ-^a-^, s = 
and the dimensionless surface charge density (Eq.(j3])) \s = gl = 0.065/a^. ^' is in units of 



^^ksTc, and L is in a units, where a is the molecular diameter (a ~ Inm for lutidine). 



large surface charge it could overcome the electrostatic repulsion. We calculated "^{L) from 
Eqs.(l8l)-(|T2l) by inserting the solutions of Eq.(lT3l) and linearized Eq.(fT^ with the boundary 
conditions (1151) . The result is shown in Fig. 3 for a particular model system. Indeed, the 
potential is repulsive far from the critical point because the electrostatic repulsion dominates, 
becomes attractive and again repulsive when the critical temperature is approached. 

The above theory is derived from the microscopic statistical mechanical description by a 
systematic coarse-graining procedure. We neglected any difference in the chemical nature of 
the cation and the anion. Coupling between the excess concentration of the solvent, "^liz), 
and the charge density, (f){z), results first from the coupling between "^ilz) and '&2{z) in 
Eq. ffTOl) . originating from the large difference in the solubilities of the hydrophilic ions in the 
two components of the solvent, and next from the coupling between '&2{z) and (f){z) of the 



entropic origin (Eq. fll2l) ). Very recently similar behavior of "^{L) was obtained in Ref. 12| 



The change of the adsorption preference in Ref. 121] results from different solubilities of the 



anion and the cation in water. Further studies are necessary to verify which mechanism 



plays the key role for the experimental results reported in Ref. ^\ 
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